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ABSTRACT 

' Atmospheric turbulence is a serious problem for ground-based interferometers. It places tight limits on both 
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sensitivity and measurement precision. Phase referencing is a method to overcome these limitations via the use 
of a bright reference star. The Palomar Testbed Interferometer was designed to use phase referencing and so 
. can provide a pair of phase-stabilized starlight beams to a second (science) beam combiner. We have used this 
' capability for several interesting studies, including very narrow angle astrometry. For close (1-arcsecond) pairs 
■"nJ" ■ of stars we are able to achieve a differential astrometric precision in the range 20-30 micro-arcseconds. 
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1. INTRODUCTION 
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-4-^ ■ Although they provide very high angular resolution, ground-based stellar interferometers suffer from very limited 
^ I sensitivity. This is a consequence of atmospheric turbulence which corrupts the incoming stellar wavefront on 
very short timescales and over small distances. Typical numbers for atmospheric coherence length and time are 
Vq = 10cm and Tq = 2ms at visible wavelengths. An interferometer must collect enough photons in a coherence 
. volume (~ rQTQ) to detect fringes. Optimistically assuming a system throughput of 10%, this would place the 
H ' S/N~3 tracking limit at ~ 8th magnitude, assuming no detector noise. 

There are several ways to overcome this sensitivity limit. The most obvious ~ and most expensive - is to 
put the interferometer in space, ^ or at least at some location with favorable atmospheric properties such as 
Antarctica.^ It is also crucial to maximize system throughput and reduce other noise sources such as detector 
read-noise.^ It should also be noted that both tq and tq scale approximately as A^^^ and thus going to longer 
wavelengths will rapidly improve the situation. This is why the current generation of interferometers often work 
in the near-IR. 

In recent years, much work has been done to increase the apparent ro via the use of adaptive optics^ (AO); 
this has been successful at both the Keck Interferometer'^ and the VLTI.^ The use of AO-corrected apertures 
provides a substantial improvement in collecting power, i.e. by a factor of ~ S" x (Z?/ro)^, where S is the AO 
Strehl ratio and D is the size of the AO-corrected aperture. For a large aperture and a typical AO Strehl one 
might expect an improvement by a factor of 100-300, as long as a natural or laser guide star is available. 

However, it is not clear that such an improvement will be sufficient for certain areas of research, e.g. for 
imaging faint extra-galactic objects or for narrow-angle astrometry.^ An additional approach is to find a way 
to increase the effective atmospheric coherence time, tq. This can be done by simultaneously observing a bright 
reference star and the science target using two separate interferometric beam-combiners, and using the reference 
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star to determine the fringe motion introduced by the atmosphere. This is known as phase referencing. By using 
a sufficiently fast control loop (~ 100 Hz) it is possible to remove much of the fringe motion in real time, and 
hence provide an artificially long coherence time for the beam combiner used to observe the science target. In 
principle one could consider extending the maximum allowed integration time from milliseconds to seconds or 
even minutes, with a corresponding increase in sensitivity. Of course, it should be noted that phase referencing, 
like AO, requires a bright reference. 

In this paper wc discuss recent progress in the area of phase referencing* as applied to very narrow-angle 
differential astrometry.^ In particular, we have used the Palomar Testbed Interferometer in a phase-referencing 
mode that allows us to achieve very high precision relative astrometry between 0.1-arcsecond pairs of stars; 
such astrometry would be useful in searching for planets orbiting one of the components in the system. This 
instrumental configuration can also be used to obtain simultaneous high spectral and spatial resolution data 
through the use of double- Fourier interferometry.^*^ 

The Palomar Testbed Interferometer^^ (PTI) is a long-baseline infrared interferometer located at Palomar 

Observatory near San Diego, California. It typically operates in the H (1.6/im) and K (2.2/im) bands. It has 
three 40-cm diameter apertures separated by 110-m and 80-m baselines, giving a maximum angular resolution 
of ~ 3 milli-arcseconds. 



2. PHASE REFERENCING 

Atmospheric turbulence above the interferometer can be characterized by a distance over which the wavefront 
remains flat. This distance is known the atmospheric coherence length. For typical observing sites, we have 
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In addition, since the atmosphere is not static, the wavefront errors change on timescales of tq ^ rg/v where 
11 is a a wind speed. The resulting phase perturbations obey an /"^^^ power law. An example of measured fringe 
motion is given in Figure 1. The fringe motion will tend to "smear" the fringe, reducing fringe visibility and 
hence signal-to-noisc ratio. 

The phase fluctuations introduced by the atmosphere are correlated over small angles, characterized by the 
angular isoplanatic angle 9o such that 
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is a parameter representing the turbulence amplitude, fc = 27r/A is the wavenumber of the light and (f) is the 
zenith angle of the observation. 

If one has available an estimate of the fringe motion induced by the atmosphere, it is sometimes possible 
to correct the phase and recover the fringe signal. If the fringe is detected using an integrating detector, or a 
detector with a non-zero read-noise, it is necessary to apply the phase correction in real-time, i.e. by rapidly 
adjusting the instrumental delay using a fast delay line or similar device. The phase reference can be obtained 
in many different ways; it can be provided by a nearby reference star (located within the isoplanatic patch, ~ 30 
arcseconds in the near-IR), or it can be the science target itself, observed at a different wavelength or with a 
different beam-combiner (which may be optimized for phase measurement). It may even be possible to provide a 
phase reference by measuring some fraction of the atmospheric pathlength changes using a laser. In any case, 
the phase referencing system needs a fast control loop to sense and correct the atmospheric pathlength errors 
on millisecond timescales. The effect of atmospheric turbulence, with and without active phase referencing is 
shown in Figure 2. 
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Figure 1. Power spectral density of tlie measured fringe pliEise, obtained using a fringe tracker at PTI, witli and witliout 
phase referencing. Tlie atmosphere induces fringe motion that has a power-law dependence {P{f) oc f~^^^) 
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Figure 2. The effect of integration time on measured fringe visibility, with and without phase referencing. Model 
parameters were appropriate for PTI, with a fringe tracker update rate of 100 Hz. On timescales shorter than the fringe- 
tracker closed-loop bandwidth, phase referencing has little effect. However, it does remove the fringe motion at low 
frequencies, which would - if left uncorrected - completely wash out the fringe. 



Clearly, phase referencing opens up a new regime for interferomctry, as beam-combination instruments can 
begin to work with exposure times measured in seconds. This will allow the use of instruments that are not 
necessarily optimized for short-exposure phase measurements but rather high spectral resolution or high-precision 
visibility measurements. This capability may also be applicable to multi-aperture combination,^'^ in that for an 
N-aperture interferometer it is only necessary to fringe-track on N baselines in order to keep the array "phased" . 
The full set of N(N-l)/2 baselines can then be combined in a different, dedicated beam-combiner. 

3. ASTROMETRY 

One application that can take advantage of phase referencing is narrow-angle astrometry. In this case one is 
interested in measuring the angle between a close pair of stars to very high precision; such astrometry can be 
used both to determine orbital parameters (which can be combined with radial velocity and photometry data 
to determine component masses, distances and luminosities) to high precision, and to search for planets in the 
system. 

PTI is unique among current interferometers in that it has a "dual-star" configuration, which makes it possible 
to simultaneously track two close (separation < 60 arcseconds) stars. To achieve this, PTI is equipped with two 
separate optical trains and beam-combiners that share the same apertures. PTI is also particularly suited to 
phase-referencing in that it is capable of fringe-tracking at bandwidths high enough to follow atmospheric motion, 
and hence one beam combiner can track and correct the atmosphere while the other can be used for high-precision 
measurements. 

For close pairs of stars, astrometric errors introduced by the atmosphere will tend to be common-mode, and 
hence will subtract out of a differential measurement. Shao & Colavita'' quantified this effect (in arcseconds of 
astrometric error) as 

C7„ = 54OB-2/30t-i/2 (4) 

where B is the baseline in meters, 6 is the angular separation (in radians) of the stars, and t the integration time 
in seconds. 
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Figure 3. The expected narrow-angle astrometric performance in arcseconds for the phase- referenced fringe-scanning 

approach, for a fixed scan rate. There arc three primary sources of astrometric error in this method: angular anisopla- 
natism,*^ temporal decoherence,*^ and photon noise. Also shown is the magnitude of the temporal decoherence effect in 
the absence of phase referencing, illustrating why stabilizing the fringe via phase referencing is necessary. 

From Equation 4 it is clear that at small separations the atmospheric error becomes negligible. This motivated 
us to modify PTI so as to allow us to study binaries with separations on the order of 0.1-1 arcseconds, i.e. larger 




Figure 4. The instrument configuration for piiase-referenced narrow-angle astrometry. One beam-combiner is used to 
traclc and correct for atmosplieric turbulence; corrections it derives are applied to both delay lines. The second beam 
combiner scans over ~ 200/im in differential optical path, tracing out the double-fringe packet produced by a close binary 
star. 

than the nominal interferometric field of view (given by A/AA x X/B), but smaller than the 30-60 arcseconds 
previously considered for narrow-angle astrometry.'' The probability of finding a background reference star 
in this reduced field is very small, but outweighed by the fact that binary stellar systems with these apparent 
separations are quite common. 

For a binary separation of 0.1 arcsecond and an interferometric baseline of 100 meters, the predicted error 

due to the atmosphere is ^ 0.2/iarcsecond in an hour. Of course, there are additional sources of astromctric 
error, including photon noise and potential systematic errors associated with the instrument. In particular, if 
the measurements of the two stars are not truly simultaneous - as in this case, see below - there will be a small 
error associated with fringe motion during the time between the measurements. Nevertheless, we estimate that 
PTI should be able to obtain an astrometric precision at the level of 10/;iarcseconds in an hour of observation 
(Figure 3). 

3.1. Instrument Configuration 

We configured PTI to observe close (separation less than 1 arcsecond) visual binaries, so-called "speckle binaries" . 
In this configuration (Figure 4), the incoming starlight was split 70/30 between a fringe tracker running at 100 
Hz and a second, astrometric beam combiner. The astrometric beam combiner did not track the fringe but rather 
it measured the intensity of the combined stellar beams while the relative optical path between the arms of the 
interferometer was modulated (using a differential delay line) in a triangle-wave pattern with an amplitude of 
~ lOO/xm and a period of ~ IHz. An example of this data is given in Figure 5. 

3.2. Data Processing 

We extract the relative astrometric positions from data such as that shown in Figure 5 as follows. First, we 
construct a search grid in differential R.A. and declination over which to search. For each point in the search 
grid we calculate the expected differential delay based on the interferometer location, baseline geometry, and 
time of observation for each scan. A model of a double-fringe packet is then calculated and compared to the 
observed scan to derive a value; this is repeated for each scan, co-adding all of the values associated with 
that point in the search grid. The final surface as a function of differential R.A. and declination is thus 
derived. The best-fit astrometric position is found at the minimum-x^ position, with uncertainties defined by 
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Figure 5. Measured intensity in the detector as a function of differential optical patli, for successive scans of tfio speckle 
binary system HD 171779. Each scan takes 1.5 seconds to acquire. The fringe tracker was locked on to the bright star 
(around 0), while the second star produces a fringe pattern which starts at -60 /im and moves due to baseline rotation. 
Although the second fringe pattern is relatively faint, the effect of coherently co-adding 500-2000 scans produces a high 
signal-to- noise ratio in the final astrometric measurement. 

the appropriate contour - which depends on the number of degrees of freedom in the problem and the value of 
the x^-minimum. Because the data was obtained with a single-baseline instrument, the resulting error contours 
are very elliptical, with aspect ratios at times > 10. 

One potential complication with fitting a fringe to the data is that there are many local minima spaced at 
multiples of the operating wavelength. If one were to fit a fringe model to each scan separately and average (or 
fit an astrometric model to) the resulting delays, one would be severely limited by this fringe ambiguity (for a 
110-m baseline interferometer operating at 2.2/im, the resulting positional ambiguity is ~ 4.6 milli-arcseconds) . 
However, by using the x^-surface approach, and co-adding the probabilities associated with all possible delays 
for each scan, the ambiguity disappears. This is due to two things, the first being that co-adding simply improves 
the signal-to-noise ratio. Second, since the observations usually last for an hour or even longer, the associated 
baseline change due to Earth rotation also has the effect of "smearing" out all but the true global minimum. 
The final x^-surface does have dips separated by ~ 4.6 milli-arcseconds from the true location, but they only 
show up at the 40-cr contour level. 

3.3. Preliminary Results 

We have obtained exploratory data on 25 binaries, and have intensively followed the binary system HD 171779 
in order to evaluate the performance and repeatability of our system (Figure 6). HD 171779 is a pair of K-giants 
in a long-period orbit, such that the apparent motion is nearly linear over the few months of observation. A 
typical observation of this pair takes an hour, and produces a formal l-a error ellipse approximately 5 x 100 /x- 
arcseconds in size; the smaller axis being aligned with the mean baseline orientation (which in turn depends on 
the exact timing of the observation). A linear two-dimensional (taking into account the full shape of the error 
ellipse) fit of 15 nights of data produces a reduced a/x^ = 4.7, indicating that the formal errors underestimate 
the true errors by a factor of 4.7. We note that fits to 1-D projections (e.g. along the declination or R.A. 
axis) show similar values of reduced x^) indicating that the error scaling is uniform. Hence we conclude that 
the night-to-night repeatability is 24 and 470 /i-arcseconds in the minor and major axes respectively. We note 
that we have recently improved aspects of the instrument metrology system and data processing, and expect to 
achieve the 10 /x-arcsecond goal in the next few months. 
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Figure 6. Measured differential astrometry between the components of the speckle binary HD 171779, together with a 
best-fit orbital motion curve. The contours represent the formal l-cr and 4-a uncertainties. Note that the single-baseline 

configuration used to obtain this data is sensitive primarily in the declination direction, and therefore the Declination- 
axis has been magnified by a factor of 10 compared to the R.A. axis. This data implies a night-to-night astrometric 
repeatability of 24 pi-arcseconds in the direction of greatest sensitivity (approximately aligned to the declination axis), 
and a somewhat less impressive 470 /it-arcsecond repeatability in the orthogonal direction. 



4. DISCUSSION 

Ground-based interferometers are severely challenged by the effects of atmospheric turbulence. The short inte- 
gration times required in the absence of any active compensation limits intcrfcrometry to all but the brightest 
astronomical sources, and limits the types of observations that can be done (spectral resolution, multi-aperture 
combination etc). One solution is to divide the problem into two parts which can be optimized separately. 
A phase-tracking beam-combiner can use cither a nearby reference source or the science target itself if it is 
sufficiently bright - to stabilize the phase of the incoming stellar wavefront. A second, science-optimized beam 
combiner can then obtain desired data at a relatively leisurely pace. 

We have used a phase-referencing approach to observe close pairs of binaries (separations less than 1 arc- 
second) in order to obtain very high precision differential astrometry. Over 15 nights of data on the binary 
HD 171779 we find a night-to-night repeatability of approximately 24 /x-arcseconds in one axis. This level of 
astrometric precision is sufficient to begin a search for Jupiter-mass planets around the components of nearby 
binary stellar systems, given that simulations"'^* indicate there are regions where planets can persist indefinitely 
(i.e. sufficiently close to one of the stellar components). A dedicated astrometric survey could quickly answer 
the question of whether or not planets form in these comparatively close systems. We have begun such a survey, 
PHASES ("Palomar High-precision AStrometric Exoplanet Search") and hope to have results in the next few 
years. 
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Figure 7. Phase space explored by radial velocity and astrornetry. The figure shows the areas in the senii-niajor axis 
vs. companion mass plane explored by a notional astrometry project ("PHASES"), for each specific target in our sample; 
astrornetry will find any planet located above and to the left of the arrows. This accounts for the actual mass and 
distance to the binaries in question. Also shown are a number of planets known to exist in wide (2 arcsec separation) 
binary systems. The precision for the RV sample is 20 m s~^ which is a realistic value given the additional photon noise 
and systematics present for compact binaries. For binary separations of 10 AU, companions in orbits inside approximately 
1 AU are expected to be stable; the stability zone grows with increasing binary separation. 
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